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PREFACE 


The  study  described  in  this  report  was  conducted  during  the  period 
July  1977  to  March  1979  for  the  U.  S.  Army  Engineer  District,  Tulsa,  by 
the  U.  S.  Army  Engineer  Waterways  Experiment  Station  under  the  general 
supervision  of  Messrs.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory, 
and  M.  B.  Boyd,  Chief  of  the  Hydraulic  Analysis  Division.  A  preliminary 
report  presenting  the  results,  conclusions,  and  recommendations  from  the 
investigation  was  presented  to  the  Tulsa  District  in  March  1979 • 

The  unsteady  flow  mathematical  model  calculations  and  tow  squat 
model  development  and  application  were  accomplished  by  Mr.  Carl  Huval. 

Mr.  Huval  also  prepared  this  report  for  publication.  Special  thanks 
should  go  to  Dr.  Billy  Johnson  and  Mr.  Paul  Senter  for  help  in  model 
application  and  data  preparation.  The  writer  appreciates  the  help  of 
Dr.  Garbis  Keulegan  and  Mr.  M.  B.  Boyd  and  Mr.  Frank  A.  Herrmann,  Jr., 
for  their  careful  review  and  suggestions  on  this  report. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  study 
and  the  preparation  and  publication  of  this  report  were  COL  John  L. 
Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customsry  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (Si)  units  as  follows: 


Multiply _  _ By _  _ To  Obtain 


cubic  feet  per  second 

0.02831685 

cubic  metres  per  second 

cubic  yards 

0.761*55^9 

cubic  centimetres 

feet 

0.301*8 

metres 

feet  per  second 

0.301*8 

metres  per  second 

feet  per  second  per  second 

0.301*8 

metres  per  second  per  second 

miles  (U.  S.  statute) 

1.60931*1* 

kilometres 

square  feet 

0.09290301* 

square  metres 

tons  (2000  lb,  mass) 

907.185 

kilograms 

i  $  nSftt 


LOCK  APPROACH  CANAL  SURGE  AND  TOW  SQUAT  AT 
LOCK  AM)  DAM  17 ,  ARKANSAS  RIVER  PROJECT 


Mathematical  Model  Investigation 


PART  I :  INTRODUCTION 

Project  Description 

1.  Lock  and  Dam  17  (also  known  as  Choteau  Lock  and  Dam)  is 
located  on  the  Verdigris  River  a  few  miles  downstream  of  the  head  of 
navigation  of  the  Arkansas  River  project  at  Catoosa,  Okla.  (Figure  1). 

The  navigation  features  of  the  project  provide  for  a  9-i“t-deep*  channel 
hy  a  system  of  locks  and  dams  on  the  Verdigris  and  Arkansas  Rivers.  The 
lock  chamber  is  110  by  600  ft  long.  Minimum  channel  width  of  150  ft  was 
dredged  for  the  Verdigris  reach  of  the  project. 

2.  After  considering  several  plans.  Dam  17  was  built  on  the 
Verdigris  River  at  mile  9-8,  while  the  lock  was  located  in  a  canal  about 
3400  ft  east  of  the  Verdigris  River  channel.  Figure  2  presents  the  lay¬ 
out  of  the  pertinent  features  of  the  project.  The  normal  upper  pool  at 
L&D  17  is  at  el  511-0**  and  the  normal  lower  pool  established  by  Dam  l6 
(Webbers  Falls)  is  at  el  490.0.  The  lock  has  a  normal  lift  of  21  ft 
and  a  maximum  lift  of  24  ft.  The  dam  and  spillway  as  built  consisted 
of  three  60-  by  27-ft  tainter  gates  with  the  spillway  crest  at  el  485.0 
and  an  earth  embankment  extending  from  the  spillway  to  the  lock.  Details 
of  the  canal  entrance  and  spillway  location  were  developed  with  the  aid 
of  a  physical  model  study  at  the  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  (Franco,  Glover,  and  Melton  1970). 

3-  The  upstream  approach  canal  to  Lock  17  is  nominally  150  ft 
wide  and  9  ft  deep  for  a  distance  of  about  5560  ft  as  shown  in  Figure  2. 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level  (msl). 
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Details  ox-  the  canal  and  river  dimensions  are  given  in  Table  1.  The 
canal  includes  two  wide  reaches  to  aid  navigation  and  reduce  surge  ef¬ 
fects  at  both  ends  of  the  150- ft  canal.  A  lll6-ft-long  reach  with  300-ft 
bottom  width  is  provided  at  the  canal  entrance  from  the  Verdigris  River. 

A  521-ft-long  reach  by  300-ft  bottom  width  was  also  excavated  near  the 
upstream  lock  approach.  The  canal  includes  two  transitions  from  the 
nominal  dimensions  of  150-ft  bottom  width  to  300-ft  bottom  width.  The 
upstream  1:2  transition  contracts  the  canal  as  tows  travel  downbound  and 
the  downstream  1:10  transition  contracts  the  canal  as  tows  travel  up- 
bound.  The  150-ft-wide  canal  conforms  to  the  authorizing  document  mini¬ 
mum  width  for  the  project.  The  design,  however,  included  provision  for 
future  widening  as  warranted  by  navigation  traffic. 

1.  iMost  of  the  towboat  traffic  operating  on  the  Verdigris  River 
is  in  the  2000-hp  class  but  regular  runs  are  being  made  by  a  few  boats 
having  up  to  1200  hp.  Tows  longer  than  600  ft  experience  difficulty 
negotiating  the  Verdigris  channel  bends.  Most  of  the  loaded  tows  on  the 
Verdigris  are  102  to  108  ft  vide  by  7-  to  8-1/2-ft  draft  and  slightly 
under  600  ft  long,  which  is  the  maximum  size  that  can  be  accommodated  in 
the  lock  without  breaking  tows.  Figure  3  snows  a  typical  canal  cross 
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Figure  3.  Tow  in  canal 


section  with  an  8-1/2-ft  draft  by  105-ft  beam  tow  in  the  canal.  There 
is  only  about  1/2-  to  1-ft  clearance,  depending  on  position  along  the 
canal,  on  each  edge  of  the  canal  as  shown  in  Table  1.  The  tow  cross 
section  is  a  major  fraction  of  the  total  canal  section. 

Navigation  Problems 


5-  Navigation  problem  reports  from  towboat  operators  have  been 
received  by  the  lockmaster  on  a  number  of  occasions.  The  large  frontal 
area  of  most  of  the  tows  in  comparison  with  the  canal  cross-sectional 
area  has  been  observed  to  cause  a  towboat  to  squat  as  much  as  1-1/2  to 
2  ft  below  its  static  floating  position,  depending  on  tow  size  and  speed. 
Information  from  the  field  indicates  that  tow  grounding  problems  occur 
mainly  while  maneuvering  through  the  canal  transitions  and  when  acceler¬ 
ating  in  the  canal.  However,  most  of  the  larger  towboats  can  ground  on 
the  bottom  going  steady  ahead  in  the  narrow  canal  reach.  Problems  with 
grounding  on  the  bottom  by  loaded  tows  occur  both  unbound  and  downbound. 
The  unbound  tows  report  problems  at  a  point  which  is  about  600  ft  down¬ 
stream  of  the  beginning  of  the  1:10  transition  (Figure  2).  Downbound 
tows  ground  at  a  location  about  200  ft  downstream  of  the  beginning  of 
the  1:2  transition. 

6.  The  squat  problem  can  be  further  aggravated  by  filling  the 
lock  chamber  while  a  downbound  tow  is  in  the  approach  canal.  If  possi¬ 
ble,  the  lock  is  filled  well  before  downbound  tows  reach  the  approach 
canal  to  minimize  the  effects  from  the  lock  surge.  As  traffic  increases, 
however,  this  operational  constraint  could  cause  considerable  delays. 
Field  tests  have  shown  that  filling  of  Lock  17  can  cause  as  much  as 
1.3  ft  of  drawdown  at  a  point  1  mile  upstream  from  the  lock.  Squat  or 
the  lock  filling  surge  can  cause  grounding  of  the  tow,  especially  when 
the  pool  elevation  is  low  and  the  depth  is  at  the  minimum  of  9  ft. 
Combinations  of  lock  surge  and  tow  squat  conditions  could  cause  grounding 
of  a  severity  that  could  lead  to  rupturing  a  barge  or  breaking  up  a  tow. 

7-  Hue  to  the  difficulties  from  lock  filling  surge  and  tow  squat, 
the  approach  canal  to  Lock  17  has  had  insufficient  depth  to  accommodate 


fully  loaded  tows.  The  U.  S.  Army  Engineer  District,  Tulsa,  thus 
requested  authority  to  widen  the  upstream  approach  channel  at  Lock  17 
from  150  to  300  ft  and  to  increase  the  depth  from  a  nominal  9  to  12  ft. 
Factual  data  were  not  sufficient  to  verify  this  plan  to  be  the  most 
economical  and  effective  method  to  relieve  the  problem.  It  was  suggested 
that  simply  widening  the  channel  from  150  to  300  ft  may  not  eliminate 
grounding,  particularly  at  times  when-  the  channel  depth  in  only  9  ft. 
Thus,  increasing  depths  from  9  to  lH  ft  may  be  much  more  cost-effective 
for  an  equivalent  amount  of  excavation.  It  was  desirable  that  the 
design  of  the  canal  remedial  plans  should  consider  the  quantities  of 
materials  involved  for  each  of  several  alternatives  to  get  a  relative 
idea  of  the  costs  for  each  plan. 

8.  WES  was  requested  to  conduct  a  numerical  model  investigation 
of  the  lock  surge  and  towboat  squat  in  the  Lock  17  approach  canal.  The 
numerical  model  study  was  conducted  to  determine  the  most  effective 
means  for  resolving  the  navigation  problem  and  to  minimize  the  towboat 
squat  and  lock  surges  from  filling  operations.  The  purpose  of  the  study 
was  to  assist  the  Tulsa  District  in  determining  the  minimum  canal  width 
and  depth  that  would  provide  adequate  navigation  conditions. 

Scope  of  Study 

9.  The  canal  surge  problem  was  modeled  mathematically  using  an 
available  unsteady  flow  model  (Johnson  197*0.  This  model  had  been  used 
to  study  several  open-channel  flow  problems  including  surge  conditions 
in  the  Illinois  Waterway  resulting  from  lock  filling  and  emptying 
operations.  The  region  modeled  for  this  study  includes  the  8500-ft-long 
canal  upstream  from  the  lock  and  a  1210-ft-long  reach  of  the  Verdigris 
River  to  a  point  just  upstream  of  the  old  Highway  69.  Figure  k  shows  a 
sketch  of  the  channel  system  that  was  modeled.  Lock  surge  field  tests 
were  conducted  in  September  1976  and  March  1978  and  these  data  were 
used  for  calibration. 

10.  Squat  effects  in  the  canal  were  computed  in  separate  calcu¬ 
lations  assuming  the  towboat  is  moving  at  a  steady  speed  ahead  in  the 
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canal.  The  steady  squat  effects  were  modeled  using  the  approach  out¬ 
lined  by  Schijf  (1953).  The  model  study  included  comparative  surge  and 
squat  calculations  of  canals  with  the  following  depths  and  widths: 


Width 

Depth 

ft 

ft 

150 

9-12-1*1-15 

200 

9-12-1*1-15 

250 

9-12-1*1-15 

300 

9-12-1*1-15 

These  cross  sections  are  shown  in  Figure  5-  The  design  tow  size  was 
chosen  as  8-1/2-ft  draft  and  105-ft  width,  which  gives  about  893-sq-ft 
cross-sectional  area. 

WIDENING  (B.W.) 

0  t  1 50'  200'  250'  300’ 


Figure  5.  Canal  cross  section 


PART  II:  MODEL  DESCRIPTIONS 


Lock  Surge  Model 

11.  The  numerical  model  used  to  determine  lock  filling  surge 
heights  along  the  lock  canal  is  a  model  called  SOCHMJ  (Stimulation  of 
Open  Channel  Hydraulics  in  Multi-Junction  Systems)  (Johnson  197*0- 
The  development  of  the  computer  program  called  SOCHMJ  has  been  docu¬ 
mented  in  the  cited  report  and  has  been  in  general  use  by  the  Math  Model¬ 
ing  Group  for  flood  routing,  tidal  computations,  and  channel  surges. 

The  equations  solved  by  the  model  are  based  on  unsteady  flow  and  are 
derived  in  books  on  open-channel  flow  (Henderson  1966).  The  resulting 
equations  of  motion  can  be  written  as: 


Continuity:  f  ♦  |  ^  }  -  0 


(1) 


Momentum: 


3V 


3V 


3h 


+  V  —  +  st  —  +  I.Y.l...  =  0 

9t  9X  8*  2.21R*'3 


(2) 


where 


h  =  water-surface  elevation  above  mean  sea  level* 

3/3t  =  rate  of  change  with  respect  to  time 
W  =  top  width  of  water  surface 
A  =  cross-sectional  flow  area 
V  =  mean  flow  velocity 

3/3x  =  rate  of  change  with  respect  to  distance 

q  =  lateral  inflow  per  unit  distance  along  the  channel  per  unit 
time 

g  =  acceleration  due  to  gravity 
n  =  Manning's  resistance  coefficient 
R  =  hydraulic  radius 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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These  equations  are  often  referred  to  as  the  equations  of  St.  Venant. 
Sir.ce  analytical  solutions  don't  exist,  numerical  techniques  such  as 
finite  differences  must  be  employed.  The  centered  explicit  finite  dif¬ 
ference  scheme  developed  by  Stoker  has  been  utilized  in  SOCHMJ. 

12.  Since  an  explicit  finite  difference  computation  scheme  is 
utilized,  the  stability  criterion 


v  +  J7 1 <  i  - 


W  /  Ax  - 


2.21R 


4/3 


must  be  satisfied  by  the  time  and  spatial  steps,  i.e..  At  and  Ax  ,  to 
ensure  a  stable  solution  of  the  difference  equations. 

13-  Data  required  for  the  operation  of  SOCHMJ  are  read  from  com¬ 
puter  cards.  The  first  data  card  contains  basic  information  such  as  the 
total  number  of  net  points  in  the  system,  the  total  number  of  junctions 
and  branches,  and  the  time-step  employed  in  the  computations.  For  the 
application  to  the  lock  approach  canal  these  values  were  99,  7,  9»  and 
3  sec,  respectively.  Figure  4  shows  the  model  network  and  Tables  2 
and  3  give  the  schematization  data.  The  three  outer  branch  ends  are 
considered  to  be  boundary  conditions  and  not  junctions.  The  second 
data  group  contains  information  about  each  of  the  nine  modeled  branches. 
Such  information  consists  of  the  numbers  of  the  first  and  last  net  points 
of  each  branch,  the  type  of  boundary  conditions  prescribed  at  the  vari¬ 
ous  boundaries,  and  the  size  of  the  spatial  step  to  be  employed  for 
each  branch.  The  third  data  group  contains  information  about  the  junc¬ 
tions  in  the  system  being  modeled  such  as  the  numbers  of  the  branches 
joining  each  junction.  The  next  major  data  group  consists  of  the  tables 

of  geometric  data.  A  table  of  top  width,  flow  area  (hydraulic 
2 /3 

radius)  ,  and  Manning's  n  (all  as  functions  of  elevation)  must  be 
input  at  each  net  point  in  the  system.  These  tables  were  developed  from 
cross-sectional  plots  as  in  Figure  3  furnished  by  the  Tulsa  District  for 
canal  and  river  cross  sections  as  shown  in  Figure  4.  The  tables  con¬ 
sisted  of  nine  elevation  data  value  entries  to  describe  the  channel  geom¬ 
etry.  The  Manning's  n  used  was  0.025  for  all  the  computations. 
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14.  The  fifth  data  group  specifies  initial  values  of  the  eleva¬ 
tion  and  discharge  at  all  grid  points  on  the  first  two  tine  lines. 

Initiax  conditions  were  set  for  zero  flow  and  5H*3-ft  elevation  through¬ 
out  the  model.  Sensitivity  tests  indicated  these  initial  conditions 
have  a  very  minor  effect  on  canal  surges  due  to  initial  steady  riverflow 
of  7000  cfs  as  compared  with  the  no-flow  condition.  The  final  major 
data  group  required  by  SOCHMJ  consists  of  the  time-dependent  boundary 
conditions  which  must  be  prescribed  at  each  open  boundary.  At  such  a 
boundary,  elevations,  discharges,  or  a  rating  curve  may  be  prescribed  as 
the  boundary  condition.  For  the  study  described  herein,  a  time-discharge 
plot  of  lock  filling  furnished  by  the  Tulsa  District  was  used  on  the 
downstream  end  of  the  model  as  shown  in  Figure  6.  Adjustments  required 
in  the  lock  filling  curve  are  discussed  later  in  the  section  on  model 
calibration.  Several  trial  computations  were  conducted  to  set  the  bound¬ 
ary  conditions  on  the  upper  ends  of  the  model.  Based  on  these  tests, 
a  zero  inflow  boundary  condition  in  the  upper  end  of  the  river  and  an 
elevation  of  511-3  at  the  spillway  were  used. 

15-  Tabular  computed  output  can  be  obtained  from  SOCHMJ  at  speci¬ 
fied  stations  along  the  modeled  reach  as  well  as  at  all  or  specific 
computational  net  points.  At  such  net  points,  output  in  the  form  of 
elevations,  velocities,  and  discharges  plus  geometric  data,  if  desired, 
is  provided.  The  output  time  interval  is  variable  and  can  be  selected 
in  the  input  data. 


16.  The  squat  of  a  ship  or  tow  in  a  canal  can  be  determined  based 
on  a  one-dimensional  model  apparently  first  formulated  in  1901  by  Thiele 
(Van  de  Kaa  1978)  but  greatly  augmented  by  Schijf  (1953).  If  the  tow  is 
moving  steady  ahead  in  a  canal,  there  is  acceleration  of  the  water  parti¬ 
cle  velocity  past  the  sides  and  bottom  of  the  tow  and,  by  Bernoulli's 
equation,  a  consequent  lowering  of  the  water  level.  A  definition  sketch 
of  the  model  is  shown  in  Figure  7.  By  the  principle  of  flow 
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MOTION  RELATIVE  TO  CANAL 


MOTION  RELATIVE  TO  TOW 


CROSS  SECTION 

Figure  7.  Squat  model  definition  sketch 

vt  (\**f 

2g  2g  Z 


where 


V  =  steady  ahead  tow  speed 

X* 

u  =  return  velocity 
z  =  vertical  tow  squat 

W  =  canal  top  width  at  the  undisturbed  water  level 
A  =  canal  area  at  the  undisturbed  water  level 
a  =  statically  submerged  tow  area 

The  assumption  is  made  in  the  above  that  the  squat  and  the  canal  width 
variation  is  small  and  the  decrease  in  canal  area  can  be  approximated  by 
a  rectangle.  The  average  channel  depth  D  =  A/W  can  be  used  to  describe 


the  water  depth  for  most  practical  canal  cross  sections  using  the 
approximation. 

IT.  The  equations  can  be  reduced  to  the  following 

u2  +  2uV 


u  =  V, 


/  a  +  Wz 
t  A  -  a  -  Wz 


The  tow  squat  for  a  given  canal,  tow  size,  and  speed  can  be  calculated 
from  the  above  by  trial  for  the  squat  model.  An  initial  estimate  of  u 
is  found  by  ignoring  z  in  the  second  equation  to  obtain 


ui  =  vt(rrr) 


where 


A 

n  =  — 
a 


The  first  estimate  of  the  return  velocity  u^  is  used  in  Equation  5  to 


obtain 


t  2n  -  1 


1  2g  (n  -  1) 


The  improved  value  z is  used  for  a  refined  estimate  of  the  return 
velocity  in  Equation  6 


and  the  squat  is  further  improved  by  using 


1  /  2 
“  2g  \u; 


2  2g  \  2  2 


The  iteration  for  calculation  of  squat  and  return  velocity  is  thus  seen 
to  take  the  following  forms 


Jk+1  2g 


n  -  1  - 


n  -  1  - 


,n  -  1  - 


The  error  tolerance  was  set  to  terminate  the  iterations  by 


lZk+l  Zkl 


An  HP65  program  was  developed  to  calculate  the  squat  and  return  velocity 
using  the  calculation  procedure  outlined  above  and  is  given  in  Table  1». 
Schijf  limiting  speed 

18.  One  of  the  most  important  results  from  the  Schijf  theory  is 
that  at  certain  tow  speeds,  critical  flow  (in  the  open-channel  flow 
sense)  will  occur  adjacent  to  the  tow.  As  a  tow  approaches  this  Schijf 
limiting  speed,  tow  squat  increases  to  the  highest  level,  waves  and 
surges  create  navigation  and  canal  bank  hazards,  and  tow  resistance  be¬ 
comes  very  high.  The  limiting  speed  is  thus  an  important  upoer  speed 


parameter  and  can  be  calculated  using  the  following  explicit  formula 
given  by  Balanin  (1977) 


(9) 


19.  When  tow  speeds  approach  V  (say  V  =  0.95  V  ),  the  itera- 

x  0  x 

tion  for  squat  calculations  presented  in  Table  h  will  not  converge.  The 
following  formula  gives  the  maximum  squat  at  : 


£  u  .  cos  /,  1\  ,  3  u  .  cos  /,  1\ 

—  =  cos  —  +  — - —  (1  -  — )  -u  cos  —  +  — - —  (1  —  — ) 
3  3  \  n/  33V  n/ 


—  + 

3 


(10) 


The  return  velocity  at  V ^  can  also  be  calculated  and  is  given  by  the 
formula 


Computations  were  made  using  these  equations  by  means  of  the  HP65  pro¬ 
gram  shown  in  Table  U. 

Tow  Squat  at  Canal  Transitions 


20.  As  a  tew  maneuvers  downbouna  through  the  lock  approach  canal, 
two  transitions  are  encountered.  The  first  is  a  contraction  from  300  to 
150  ft  wide  and  farther  downstream  an  expansion  from  150  to  300  ft  wide. 
The  canal  narrows  150  ft  in  a  longitudinal  distance  of  300  ft  at  the 
first  transition  (1:2).  Several  cases  have  been  reported  of  tows  ground¬ 
ing  on  the  canal  bottom  while  maneuvering  through  this  transition. 


Downbound  tows  have  little  difficulty  in  the  canal  expansion  farther 
downstream.  On  the  other  hand,  upbound  tows  have  reported  incidents  of 
tows  "setting  on  the  bottom"  near  the  upstree~  end  of  the  1:10  transition 
(150- ft  expansion  in  a  longitudinal  distance  of  1500  ft).  The  following 
paragraphs  present  the  methodology  that  was  used  for  analyzing  the 
"supersquat"  phenomenon  when  tows  are  driving  through  canal  transitions. 
The  method  of  analysis  is  based  on  previous  work  by  Hooft  (19^9)  and 
Marchal  (1977)  as  well  as  the  review  report  by  Kolkman  (1978). 

21.  As  tows  proceed  in  the  300-ft-wide  canal  reaches  toward  either 

canal  transition,  the  speed  will  likely  approach  the  Schijf  limiting 
velocity  .  A  value  of  90  percent  of  was  the  maximum  self- 

propelled  tow  attainable  speed  achieved  during  recent  studies  by  the 
Delft  Hydraulics  Laboratory  (Van  de  Kaa  1978).  Tows  usually  proceed  at 
much  lower  speeds,  but  90  percent  was  used  for  computations.  Computa¬ 
tions  to  be  presented  later  show  that  0-9^  in  the  present  wide  (300  ft) 
canal  will  be  higher  than  the  limiting  velocity  in  the  narrow  (150  ft 
wide)  canal  reach.  The  squat  on  the  tow  will  thus  increase  in  the  nar¬ 
rowed  canal  by  a  substantial  amount,  and  the  velocity  will  exceed  the 
maximum  or  limiting  velocity  in  the  narrowed  canal.  Experiments  and 
theory  both  indicate  considerable  surge  action  at  and  near  the  limiting 
velocity.  Additional  flow  cannot  pass  the  tow  in  the  narrowed  canal 
section  and  the  resulting  flow  "choking"  will  cause  a  positive  surge  to 

be  generated  upstream  of  the  tow.  This  surge  height  will  be  a  function 
of  the  difference  between  the  tow  speed  in  the  wide  canal  and  the  limit¬ 
ing  velocity  in  the  narrow  canal  reach.  Figure  8  helps  to  clarify  the 
phenomenon  and  gives  the  definition  of  terms. 

22.  Using  the  flow  superposition  principle  as  previously  given 
for  the  steady  ahead  squat  (paragraph  l6),  the  tow  is  made  stationary 

in  the  transition  with  an  approach  velocity  V  allowing  a  quasi-steady 
flow  analysis.  The  velocity  of  the  tow  is  given  the  magnitude  of  90  per¬ 
cent  of  the  limiting  velocity  in  the  wide  canal  reach  © 
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MOTION  RELATIVE  TO  TOW 


Figure  8.  Tow  in  canal  transition 


If  this  velocity  is  greater  than  the  limiting  velocity  in  the  narrow 
canal  (2),  flow  past  the  tow  will  he  controlled  by  the  limiting  velocity 
and  a  positive  surge  of  amplitude  n  will  be  generated  at  the  bow 


© 


of  the  tow.  If 
becomes  zero 


(V  )  is  eaual  to  the  limiting  velocity  at  (2),  n 

tyQ 


■  -w. 


is  less  than  the  limiting  velocity  at  (2), 


then  the  normal  sublimit  squat  relation  applies.  As  the  tow  encounters 
the  increased  resistance  from  the  positive  surge  n  ,  the  tow  speed  will 


oecome 


less  the  (V.  )  ,  and  approach  the  value 

w© 


W©  *  °-9W 


23. 

the  given 


The  magnitude  of  the  maximum  squat  z  can  be  calculated  for 
(V  \  using  the  steady-speed  squat  model  previously 


from  the  squat  model  and  the 


described. 


Using  the  values  of 


tow  speed  in  the  wide  canal,  the  positive  surge 
using  the  surge  continuity  equation 


n 


can  be  calculated 


(12) 


where  A  and  VJ 


are  taken  at 


n  W  V  gp 
A  +  Vsn 


D  +  n 


@ 


If  n  is  a  small  quantity  with  respect  to  D  ,  the  equation  can  be 
simplified  to 


_  _D _ 

n"  VS 

2h.  It  is  evident  from  the  discussion  presented  that  the  transient 
tow  supersquat  at  the  1:2  transition  near  the  upstream  end  of  the  canal 
is  greater  than  the  steady  ahead  squat.  The  bow  of  the  tow  will  usually 
react  hydrostatically  to  the  sum  of  +  n  (defined  as  Ah),  thus 
causing  the  tow  stern  to  sink  by  that  amount  relative  to  the  undisturbed 
water  level.  This  Ah  is  the  supersquat  at  the  stern  of  the  tow. 

2p.  Tows  grounding  and  other  navigational  problems  have  also  oc¬ 
curred  at  the  1:10  transition  farther  downstream.  It  appears,  therefore, 
that  increasing  the  length  of  transition  may  not  eliminate  the  problem. 
Due  to  whe  600-ft  length  of  many  tows  using  the  canal,  it  may  be  expected 
that  superlimiting  velocity  will  occur  at  the  bow  of  the  design  tow 
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proceeding  in  the  1:10  transition.  In  general,  it  would  be  necessary 
to  make  the  canal  transitions  very  gradual  so  that  area  differences  be¬ 
tween  the  bow  and  stern  of  the  tow  cannot  lead  to  superlimiting  veloc¬ 
ities  at  the  bow.  The  present  1:10  transition,  for  example,  leads  to 
an  area  difference  of  about  600  sq  ft  between  tow  bow  to  stern.  This 
is  about  one-third  the  present  150-ft-wide  canal  cross-sectional  area 
indicating  that  the  transition  section  will  cause  a  large  impact  on  squat. 
Transition  lengths  of  three  or  more  times  the  present  Length  of  1500  ft 
would  be  necessary  to  decrease  area  changes  below  10  percent  of  the  pres¬ 
ent  canal  area. 

26.  The  analysis  presented  assumes  the  tow  is  proceeding  at  con¬ 
stant  speed  in  the  300-ft-wide  canal  reach  without  any  change  in  engine 
setting.  It  is  probable  that  some  pilots  would  attempt  to  regain  speed 
through  the  transition  by  increasing  engine  rpm.  This  would  further 
complicate  the  problem  at  a  transition  and  would  lead  to  an  increased 
tendency  for  tow  grounding. 
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PART  III:  LOCK  SURGE  MODEL  APPLICATIONS 


Calibration 


27.  Lock  surge  field  data  were  obtained  in  July  and  September 
1976  and  surge  tests  were  repeated  in  March  1978.  The  July  1976  tests 
were  made  at  a  riverflow  of  about  2L,000  cfs  and  the  September  1976 
tests  at  7200  cfs.  Four  staff  gages  were  located  at  l/2-mile  intervals 
in  the  lock  approach  canal  with  two  more  on  the  upstream  river  channel. 
Location  of  these  gages  is  shown  in  Figure  gage  1  was  located  in  the 
lock  chamber.  The  1976  tests  involved  a  series  of  three  lock  filling 
and  emptying  operations.  The  test  operations  were  conducted  so  as  to 
maximize  the  canal  surges.  Measurements  were  made  at  1-min  intervals 
for  about  an  hour.  Only  the  negative  (lock  filling)  surges  were  of 
interest  in  this  study,  since  the  problem  was  one  of  lack  of  adequate 
depth  for  navigation. 

28.  The  1978  test  was  cond  icted  at  7600-cfs  riverflow  and  con¬ 
sisted  of  one  lock  filling  operation  with  data  obtained  at  1-min  inter¬ 
vals  for  1-1/2  hr.  The  upper  pool  elevation  was  511.30  at  the  beginning 
of  the  test.  Plots  of  the  three  data  sets  indicated  only  minor  experi¬ 
mental  differences  between  the  three  test  series  during  the  first  15  min 
of  testing.  The  September  1976  and  March  1978  data  were  used  for  model 
calibration. 

29.  Figures  9j  10,  and  11  present  water-surface  elevation  time- 
history  field  data  plotted  at  three  of  the  four  canal  staff  gage  loca¬ 
tions  (gages  3,  and  5,  respectively).  The  numerical  model  computa¬ 
tions  using  the  lock  filling  curve  originally  furnished  by  the  Tulsa 
District  (Figure  6)  showed  a  less  severe  negative  surge  than  the  field 
data.  The  maximum  difference  at  gage  3,  for  example ,  was  about  0.2  ft. 
When  the  lock  filling  curve  in  Figure  6  was  integrated  for  the  6-min 
lock  filling  time  and  compared  with  the  known  volume  of  the  lock 

(21  x  110  x  670  ft),  a  17.8^  percent  difference  was  noted.  As  a  result, 
the  original  lock  filling  curve  was  adjusted  by  increasing  all  ordinates 
by  17. 8U  percent.  The  revised  model  computations  after  lock  filling 
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Figure  11.  Model  calibration,  gage  5 
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rate  adjustment  are  also  shown  in  Figures  9,  10,  and  11,  and  they  indi¬ 
cate  much  better  agreement  with  the  field  data.  Differences  between 
field  data  and  computations  at  gage  5  after  about  11-1/2  min  are  attrib¬ 
uted  to  nonrealistic  spillway  boundary  conditions.  The  primary  canal 
surge  problem  centers  around  the  maximum,  initial  negative  surge.  For 
the  purpose  of  this  study,  only  computation  results  during  the  initial 
wave  propagation  through  the  canal  for  times  less  than  11-1/2  min  were 
used  for  further  analysis. 

Existing  Conditions 


I  i 


Surge  in  existing  canal 

30.  Computation  results  for  the  existing  canal  can  be  used  to 
ascertain  the  extent  of  the  lock  filling  surge  problem  in  the  present 
canal.  Several  plots  are  presented  and  discussed  in  the  following 
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Figure  12.  Surge  profiles  without  tow 

paragraphs  to  portray  the  present  problem  and  lay  the  basis  for  analyzing 
the  effects  of  canal  modification  alternative  plans.  Figure  12  presents 
surge  profiles  along  the  approach  canal  at  1,  3,  6,  and  9  min  after 
sta:t  of  lock  filling.  The  negative  surge  front  tends  to  flatten  daring 
propagation  upstream.  The  tail  of  the  front  can  be  quite  complex,  due 
to  the  variation  of  lock  filling  rate  and  surge  interaction  with  the 
complex  geometry  in  about  the  first  1500  ft  upstream  from  the  lock 
(branches  8  and  9)-  The  surge  magnitude  tends  to  decrease  during  up¬ 
stream  propagation.  The  speed  of  surge  propagation  is  about  l6  ft/sec 
which  is  close  to  the  theoretical  celerity  of  VgD  . 

31.  An  important  and  useful  method  of  portraying  canal  surge  is 
by  plotting  the  maximum  negative  surge  profile  along  the  canal  as  pre¬ 
sented  in  Figure  13.  The  minimum  elevation  of  this  profile  indicates 
the  location  of  the  maximum  surge  (except  in  the  immediate  vicinity  of 
the  lock)  about  2900  ft  upstream  of  the  lock  gates.  This  point  is  in 
the  upstream  end  of  the  1:10  transition,  about  200  ft  downstream  from 
the  smallest  150-ft-wide  canal  cross  section.  The  maximum  surge  is 
about  l.l»5  ft.  A  tow  on  standby  near  that  position  would  bump  or 
strike  the  bottom  as  the  surge  passes. 
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riSTANCe  UPSTREAM  of  LOCK  IN  1000  FT 


Figure  13.  Maximum  surge  profile  without  tow 
Effects  of  tow  in  canal 

32.  In  some  circumstances,  lock  filling  operations  may  he  neces¬ 
sary  with  a  tow  standing  hy  in  the  canal.  Computat:  ns  were  made  with  a 
105-ft-wide  hy  600-ft-long  tow  centrally  located  at  a^out  6700  ft  up¬ 
stream  from  the  lock.  The  assumption  was  made  that  the  downhound  tow 
pilot  would  choose  to  stand  hy  inside  the  narrow  150-ft-wide  canal,  hut 
as  far  from  the  lock  surge  as  possible,  while  waiting  for  the  lock  to 
fill.  The  tow  was  assumed  to  respond  to  the  surge  wave  hy  following  the 
free  surface  in  a  quasi-static  manner.  This  is  consistent  with  the  long 
wave  approximation  which  is  the  basis  of  the  lock  surge  model.  A  small 
space  step  is  necessary  to  resolve  the  tow  effects  on  the  numerical 
surge  model.  A  total  of  nine  space  steps  were  used,  that  is,  there  was 
a  different  A.<  in  each  of  the  nine  branches  (Table  2).  The  smallest 
Ax  was  in  tne  branch  simulating  the  tow. 

33.  The  method  of  introduction  of  the  tow  effects  on  the  negative 
surge  was  based  on  the  "sunken  tow"  concept.  The  tow  is  assumed  to  act 
essentially  as  an  abrupt  contraction  at  the  bow  and  an  abrupt  expansion 


at  the  stern.  Qualitatively,  this  is  equivalent  to  decreasing  the  canal 
cross-sectional  area  by  the  tow  cross-sectional  area.  Further  elabora¬ 
tion  of  this  concept  is  given  in  the  paper  by  Kolkman  (1978).  For  appli 
cation  to  the  canal  surge  computations,  the  tow  cross-sectional  area 
was  subtracted  from  the  canal  area  at  the  bottom  of  the  canal. 

3l4.  Figures  lU  and  15  present  two  time-histories  of  surge  passage 


TIME  IN  MINUTES 


Figure  lU.  Effect  of  tow  in  canal  at  gage  U 

downstream  of  the  tow  bow  and  at  the  midtow  cross  section.  Two  tow 
drafts  of  7-1/2  and  8-1/2  ft  were  tried  but  apparently  had  only  minor 
effect  on  the  surge.  The  effect  of  an  8-1/2- ft-draft  tow  is  to  increase 
the  negative  lock  surge  at  and  between  the  tow  and  the  lock.  The  maxi¬ 
mum  effect  at  the  rnidto'*  location  is  about  0.2  ft. 

35.  The  effect  of  the  tow  is  more  clearly  shown  in  Figure  l6 
which  portrays  profiles  of  maximum  negative  surges.  The  maximum  draw¬ 
down  of  the  surge  occurs  at  the  bow  of  the  tow,  causing  an  increase  in 
surge  amplitude  from  1.08  to  l.Uo  ft.  The  effect  attributable  to  the 
tow  is  thus  equal  to  about  30  percent  of  the  maximum  surge  without  the 
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tow.  This  implies  that  a  tow  located  at  the  lock  surge  minimum  near 
the  beginning  of  the  1:10  transition  may  be  amplified  to  about  1.9  ft. 

36.  It  had  initially  been  planned  to  conduct  additional  computa¬ 
tions  to  study  the  relative  tow  effect  at  several  other  canal  sections. 
However,  it  was  decided  to  concentrate  more  effort  on  tow  squat  effects 
since  this  appeared  to  be  more  critical  to  the  navigation  problems. 

Effects  of  Canal  Modifications 


37-  A  total  of  15  canal  configurations  were  tested  with  identical 
boundary  and  initial  conditions.  In  all  cases,  the  initial  conditions 
used  were  el  511-3  (the  observed  elevation  during  the  field  tests) 
throughout  the  model  and  zero  discharge.  The  15  tests  represented  al¬ 
most  all  the  canal  configuration  combinations  needed  to  obtain  a  com¬ 
plete  picture  of  the  impact  of  canal  geometry  changes  (widening  and/or 
deepening)  on  the  lock  surges.  One  combination  for  the  lU-ft-deep  by 
300-ft-wide  canal  was  inadvertently  missed  during  generation  of  the 
canal  geometry  data  tables. 

38.  Figure  5  shows  the  method  of  numerically  "dredging"  the  canal 
at  a  typical  cross  section.  Note  that  the  existing  canal  is  considerably 
larger  than  the  trapezoidal  designed  canal  of  9  ft  deep,  150-ft  bottom 
width,  and  l-on-3  side  slope.  With  a  pool  water-surface  elevation  of 
511.3,  the  design  canal  at  9-  and  12- ft  depths  (relative  to  normal  pool 
el  511-0)  would  have  the  following  numerical  properties  as  compared 
with  the  actual  canal  measurements. 


Depth  Depth  at  Test  Pool 
ft  el  511-3 
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39-  A  summary  of  pertinent  geometrical  canal  data  for  the  test 
conditions  is  presented  in  Table  5*  The  canal  area  increase  due  to 
widening  and/or  deepening  is  also  tabulated  as  an  index  of  probable 
relative  volume  of  canal  dredging.  For  example,  doubling  the  present 
canal  bottom  width  to  300  ft  will  require  a  69  percent  increase  in  canal 
area. 

!»0.  The  effect  of  canal  widening  alone  is  presented  in  Figure  17 
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Figure  17.  Surge  profiles  at  end  of  lock  filling  (6  min), 

canal  widened 

showing  surge  profiles  at  the  end  of  lock  filling  (6  min).  Doubling 
the  present  nominal  150-ft-vide  bottom  width  canal  to  300  ft  would  de¬ 
crease  the  surge  amplitude  by  about  0.28  ft.  The  surge  amplitude  is  de¬ 
creased  not  only  in  the  widened  canal  reach  but  also  downstream  in  the 
reach  from  the  1:10  transition  to  the  lock.  The  surge  propagation  speed 
is  not  changed  materially  by  canal  widening. 

1*1.  The  effect  of  canal  deepening  alone  indicates  much  more  com¬ 
plex  surge  profiles  as  shown  in  Figure  18.  The  surge  propagation  speed 
is  increased  by  depth  increases  but  the  maximum  amplitude  is  only 
slightly  decreased.  The  surge  is  decreased  only  by  about  0.09  ft  by 
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Figure  18.  Surge  profiles  at  end  of  lock  filling  (6  min), 

canal  deepened 


increasing  the  present  canal  to  15  ft.  It  appears,  however,  that  deepen¬ 
ing  to  12  ft  would  actually  increase  the  maximum  surge  by  a  very  slight 
amount.  It  should  be  remembered,  however,  that  the  canal  depth  of 
12  ft  would  result  in  an  additional  3  ft  of  water  under  the  tow.  Thus, 
a  negative  surge  of  about  the  same  magnitude  as  presently  occurs  would 
not  result  in  grounding  of  a  standby  tow.  Figure  19  presents  the  effects 
of  doubling  the  present  canal  bottom  width  and  deepening  to  12  and  15  ft. 
Here  again,  the  surge  profile  demonstrates  the  minor  effect  of  deepen¬ 
ing  on  surge  amplitude.  An  increase  in  depth  from  12  to  15  ft  (25  per¬ 
cent  increase)  decreases  the  surge  by  about  0.1  ft. 

^2.  Profiles  of  minimum  canal  elevations  (or  maximum  surge  am¬ 
plitudes)  for  widened  canals  are  presented  in  Figure  20.  The  largest 
surge  amplitude  in  the  profile  (except  in  the  immediate  vicinity  of  the 
lock)  occurs  near  the  upstream  end  of  the  1:10  transition  (branch  8) 
for  all  three  canals  widened  to  less  than  300  ft.  These  transition  ef¬ 
fects  disappear  at  the  widest  canal  width  when  the  300-ft  width  is  con¬ 
stant  along  the  enuire  length  of  the  canal.  Canal  widening  effects  are 
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most  pronounced  in  the  reach  where  the  canal  is  widened.  The  surge  am¬ 
plitude  decreases  as  the  width  increases  for  most  of  the  canal  length. 
1*3.  The  effects  of  increasing  canal  depth  are  shown  in  Figure  21 


Figure  21.  Maximum  surge  profiles;  effects  of 
canal  deepening 


Surge  amplitudes  are  decreased  by  increasing  depths  in  the  downstream 
end  of  the  canal  frcm  the  1:10  transition  to  the  lock.  Surge  amplitudes 
are  increased  in  the  upstream  half  of  the  canal,  but  only  by  about 
0.2  ft.  Canal  depth  increases  are  less  effective  than  width  increases 
in  surge  reduction,  but  the  additional  depth  would  increase  the  clear¬ 
ance  under  a  standby  tew  and  prevent  grounding. 

UU.  Computational  results  with  the  canal  widened  to  300  ft  and 
at  various  depths  are  shown  in  Figure  22.  Surges  are  further  reduced 
in  the  downstream  end  of  the  canal  by  deepening  the  300-ft-wide  canal. 
The  largest  canal  cross  section,  15  by  300  ft,  provides  the  largest 
reduction  in  surge  amplitude  over  most  of  +he  canal  length.  However, 
the  surge  amplitude  for  the  largest  canal  excavated  is  still  nearly 
1  ft  in  magnitude  (0.90  ft).  This  represents  a  decrease  in  surge  ampli¬ 
tude  of  about  35  percent  frcm  the  existing  canal  surge  of  1.1*5  ft. 
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Figure  22.  Maximum  surge  profiles;  effect  of  deepening 
300-ft-wi de  canal 

Green's  Lav  and  Its  Application 

45-  It  is  possible  to  study  the  lock  filling  surge  in  the  ap¬ 
proach  canal  using  a  relatively  simple  model  of  the  phenomena.  The 
intent  is  to  find  the  relative  impact  of  channel  dimension  modifications 
on  lock  surge  and  tow  squat.  The  following  paragraphs  present  the 
rationale  and  the  theory  of  the  method. 

46.  The  bed  slope  of  the  approach  canal,  for  ail  practical  pur¬ 
poses,  is  zero.  Resist;_..^e  losses  in  the  canal  will  be  small  during 
the  propagation  of  the  initial  surge  upstream  from  the  lock  chamber  to 
the  river.  The  detail  derivation  is  presented  in  standard  references 
on  open-channel  flow  or  wave  motion  (for  example,  Keulegan  1950). 

4? •  As  Figure  23  shows,  the  negative  surge  propagates  upstream 
into  the  canal  which  has  zero  initial  velocity.  The  amplitude  of  the 
negative  surge  is  a  function  of  the  rate  of  lock  filling  and  geometrical 
channel  parameters.  The  Engineer  Manual  on  lock  surges  (OCE  1949)  sug¬ 
gests  increasing  the  channel  width  and/or  depth  to  reduce  the  surge 
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Figure  23.  Lock  surge  in  variable  canal 
sections 

height.  Guidance  on  the  relative  benefits  of  channel  modifications  is 
not  provided.  It  is  possible  to  investigate  the  relative  impact  of 
canal  geometry  by  means  of  the  well-known  Green's  law.  The  law  is  a 
natural  consequence  of  assuming  equal  surge  energy  for  the  alternative 
channel  dimensions.  The  law  states  that  the  surge  amplitude  is  propor¬ 
tional  to  the  inverse  of  both  the  square  root  of  canal  width  and  the 
fourth  root  of  the  mean  depth.  In  symbols,  the  relation  is  given  as 


a 

n 


(lit) 


where 


W, 


W 


surge  amplitude  for  any  given  variable  canal  geometry 

surge  amplitude  for  the  present  unmodified  canal 

present  canal  top  width 

variable  geometry  canal  top  width 

A^/W^  Presen^  canal  mean  depth 

variable  geometry  canal  mean  depth 


U8.  The  relative  amplitudes  can  be  normalized  to  obtain  relative 
surge  reductions  by  calculating  (a.^  -  a^/a.^  for  ease  of  visualization 
of  relative  area  increases.  Table  5  presents  values  of  relative  surge 
amplitude  reductions  using  Green's  law.  These  are  compared  with  the 
mathematical  model  results  in  Figure  2h  showing  reasonably  accurate 
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COMPUTED  RELATIVE  SURGE  REDUCTIONS.  PERCENT 

Figure  24.  Comparison  of  model  computations  with 
Green' s  law 


estimates  of  the  impact  of  canal  geometry  changes  on  the  surge  by  using 
Green's  law.  Figure  25  presents  a  generalized  plot  using  Green's  law  to 
portray  the  surge  amplitude  reduction  from  canal  modifications.  The 
plot  shows  that  canal  widening  reduces  surge  amplitude  more  than  an 
equivalent  area  increase  by  canal  deepening. 


as 


/7A 


y/r 


LEGEND 


EQUI  DEPTH 


I  I'S 


relative  area  increase,  percent 

Figure  25.  Relative  lock  surge  reduction  based  on 

Green ' s  law 

J*9.  Green's  law  assume;,  equal  surge  energy  for  canal  geometry  al¬ 
ternatives.  This  is  a  reasonable  assumption  in  gradually  varied  canal 
geometries  such  as  in  the  1:10  transition.  When  canal  area  changes  are 
not  gradual,  the  surge  analysis  for  abrupt  expansions  and  contractions 
should  be  used  as  presented  by  Keulegan  (1950)  and  in  the  Engineer  Manual 
on  canal  surges  (OCE  1 9^9).  This  will  result  in  higher  surges  than  those 
determined  by  Green's  law  for  equivalent  canal  area  changes. 


Summary  of  Lock  Surge  Results 


50.  Maximum  negative  surge  amplitudes  were  abstracted  from  the 
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15  computational  runs  and  are  given  in  Table  6.  Data  are  shown  at  the 
end  of  lock  filling  (6  min)  and  the  maximum  values  from  surge  profiles 
are  given.  The  maximum  surge  magnitudes  near  the  upstream  end  of  the 
1:10  canal  transition  are  shown  plotted  in  Figure  26  as  a  function  of 


Figure  26.  Summary  of  results  of 
model  computations 


canal  area.  As  expected,  surge  amplitude  is  decreased  by  increasing  the 
canal  area.  Doubling  the  present  canal  cross-sectional  area  (bottom 
width  of  about  350  ft)  without  deepening  will  decrease  the  surge  to  about 
1  ft .  The  rate  of  decrease  of  surge  is  greater  due  to  canal  width  in¬ 
creases  than  by  canal  deepening.  Increases  of  only  canal  width,  however, 
do  not  increase  the  bottom  clearance  under  the  tow  which  is  more  critical 
at  the  time  of  maximum  surge. 


40 


PART  IV:  TOW  SQUAT  APPLICATIONS 


Steady  Ahead 

51.  Computations  were  made  using  the  Schijf  squat  model  as 
described  in  PART  II  and  are  summarized  in  Table  7-  Plotted  squat  re¬ 
sults  of  these  calculations  will  be  discussed  in  subsequent  paragraphs. 

All  computations  were  made  using  the  8-1/2- ft  draft  by  105-ft-wide  design 
tow. 

52.  Tow  squat  increases  as  the  square  of  the  tow  speed  as  shown 

in  Figure  27.  This  plot  also  presents  the  effect  of  widening  the  present 
canal  and  the  Gchijf  limiting  speed  and  the  90  percent  limiting  speed 
for  each  canal  configuration.  Squat  is  greatly  increased  when  tow 
speeds  become  larger  than  90  percent  of  the  limiting  speed.  Laboratory 
and  field  tests  have  indicated  that  self-propelled  tows  cannot  exceed 
and  are  usually  operated  from  about  50  to  90  percent  of  .  The 
squat  in  the  present  canal  at  0.9V^  is  about  0.7  ft  which  allows  only 
about  0.1-ft  clearance  underneath  the  tow  for  a  pool  elevation  of  511-3- 
Widening  the  canal  increases  the  limiting  tow  speed.  If  it  is  assumed 
that  tows  will  proceed  at  the  highest  possible  speed  of  about  90  percent 
of  ,  then  increasing  the  canal  width  will  not  alleviate  the  tow 
grounding  problems  caused  by  squat. 

53-  The  effects  of  increasing  the  canal  depth  are  shown  in  Fig¬ 
ure  28.  Squat  at  limiting  tow  speed  increases  more  rapidly  with  canal 
depth  increases  than  with  widening.  The  12- ft  canal  depth  (below  511  msl) 
will  provide  adequate  depth  to  accommodate  the  8.5-ft  tow  draft  plus  the 
C.y-it  aquae  at  O.y'v'^  . 

5l4 .  The  effect  of  canal  width  on  the  12  '’t-deep  canal  is  presented 
in  Figure  29.  Tow  squat  is  not  a  strong  function  of  canal  width  at  90 
percent  of  limiting  speed.  The  maximum  tow  squat  is  about  1.0  ft  for  the 
12-  by  300-ft  canal  compared  with  0.9-ft  squat  for  the  12-  by  150-ft 
canal.  Either  would  give  adequate  clearance  underneath  the  tow  bottom. 

55-  Increasing  canal  area  will  allow  the  tows  to  proceed  at  a 
higher  speed  than  is  now  possible  in  the  existing  canal  which  will  be 


TOW  SPEED,  PT/SCC 


1  ? 


7if- 

'•  r 


TOW  SPEED,  FT/SEC 

Figure  29.  Effect  of  canal  width  increase  on 
tow  squat  at  12- ft  depth 

of  benefit  to  the  towboat  operations  since  this  will  reduce  transit 
time  through  the  constricted  canal  reach.  Relative  possible  tow  speed 
increase  is  shown  plotted  in  Figure  30  as  a  function  of  relative  canal 
area  increase  and  is  presented  in  Table  7-  The  graph  shows  that  increas¬ 
ing  canal  depth  has  a  greater  impact  on  tow  speed  than  increasing  canal 
width. 

56.  Squat  will  increase  if  tows  are  operated  at  the  same  percent¬ 
age  of  the  limiting  velocity.  The  results  of  computations  to  determine 
the  relative  squat  at  90  percent  are  shown  graphically  in  Figure  31 

and  are  presented  in  Table  7.  The  plot  shows  that  the  relative  squat 
increases  more  rapidly  by  deepening  the  canal  than  by  widening  it.  Ta¬ 
ble  7  shows,  however,  that  the  increase  in  squat  is  small  and  less  than 
the  increase  in  canal  depth;  thus,  it  would  be  more  advantageous  to 
deepen  the  canal  than  to  widen  it,  for  a  given  increase  in  cross-sectional 
area. 

57-  Limiting  the  tow  speed  in  the  canal  is  one  possible  way  of 
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Figure  30.  Effect  of  canal  size  on  relative  tow  speed  increase 

greatly  decreasing  the  risk  of  tow  grounding.  This  is  illustrated  in 
Figure  32  which  presents  the  relative  decrease  of  squat  due  to  canal 
area  increases  if  the  tow  speed  is  held  to  the  present  maximum  tow  speed 
of  0.9V£  (this  is  approximately  2-1/2  mph).  These  data  are  also  pre¬ 
sented  in  Table  7.  Canal  deepening  and  widening  have  equal  effect  on 
squat  decrease  for  a  given  increase  in  cross-sectional  area. 

Supersquat  at  Canal  Transitions 

58.  Supersquat  computations  in  the  1:2  transition  were  made  using 
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Figure  32.  Effect  of  canal  size  on  tow  squat  at  2-1/2  mph 
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present  canal  depth  is  increased  to  12  ft  without  increasing  the  width; 
there  would  be  about  a  1+0  percent  increase  at  15-ft  depth.  Increasing 
the  width  to  more  than  200  ft  is  required  to  have  a  negligible  super¬ 
squat  change  compared  with  the  present  canal.  A  canal  width  of  greater 
than  about  250  ft  is  necessary  to  reduce  the  present  supersquat  by  more 
than  one-fourth.  The  full  width  of  300  ft  without  a  canal  transition 
reduces  the  supersquat  to  essentially  the  same  magnitude  as  the  steady 
ahead  squat  of  slightly  less  than  1  ft.  Width  increases  are  more  effec¬ 
tive  in  reducing  relative  supersquat  effects  compared  with  depth  in¬ 
creases  for  equivalent  canal  area  increases. 
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PART  V:  DISCUSSION  07  RESULTS 


59-  The  lock  surge  results  indicate  that  surge  magnitudes  will  re¬ 
main  appreciable  in  any  of  the  canal  configurations  modeled  compared 
with  the  existing  canal.  Surge  magnitudes  can  only  be  reduced  from 
about  1.5  ft  to  slightly  less  than  1  ft.  A  surge  basin  (or  other 
possibly  expensive  remedial  works)  would  be  necessary  to  reduce  the  lock 
surge  to  a  negligible  amount.  Increasing  the  lock  filling  time  from  the 
present  6  min  would  decrease  the  surge  amplitude,  but  this  might  cause 
unacceptable  delajrs  as  traffic  increases  in  the  future.  The  design  and 
operation  of  canal  and  lock  should  accommodate  the  expected  surge,  since 
further  surge  reduction  would  probably  be  cost-prohibitive. 

60.  Steady  squat  in  the  canal  depends  on  the  square  of  the  tow 
speed  going  through  the  canal.  While  the  grounding  problems  could  be 
eliminated  by  imposing  speed  limits,  this  is  not  usually  found  to  be 
practicable.  The  analysis  presented  in  this  report  was  made  by  assuming 
that  maximum  possible  tow  speeds  would  be  approached  by  some  pilots. 
Minimum  acceptable  canal  dimensions  were  defined  as  the  dimensions  which 
would  accommodate  the  design  tow  under  those  conditions.  Using  these 
criteria,  it  was  determined  that  the  present  canal  is  inadequate  to 
handle  the  design  tow.  The  steady  squat  computations  indicate  that  canal 
depths  of  12  ft  or  more  would  provide  adequate  bottom  clearance.  Canal 
deepening  is  more  effective  in  increasing  limiting  tow  speed  than  canal 
widening.  A  trial  canal  depth  was  thus  set  at  12  ft,  with  the  possibil¬ 
ity  of  increasing  or  decreasing  the  design  depth  above  or  below  this 
value  if  -warranted  by  further  considerations. 

61.  Super squat  computations  at  the  canal  transition  with  several 
combinations  of  canal  depths  and  widths  sr.uved  that  the  canal  width 
would  have  to  be  increased  to  at  least  200  ft  to  not  increase  present 
supersquat  effects.  A  uniform  canal  width  of  300  ft  would  eliminate  all 
supersquat  effects  at  both  canal  transitions.  The  impact  of  the  various 
alternative  canal  sizes  on  the  net  clearance  between  the  canal  and  tow 
bottom  is  the  main  factor  to  consider  for  canal  design.  It  is  necessary, 
however,  to  consider  the  possibility  of  coincidental  effects  from  lock 


surges  and  supersquat  from  tows  maneuvering  through  the  canal  transi¬ 
tions.  The  timing  of  lock  surges  and  possible  tow  movements  through 
these  transitions  suggests  that  surges  and  supersquat  co*nd  coincide, 
especially  in  the  upstream  transition.  An  allowance  of  about  1.0  ft  due 
to  surge  effects  at  and  in  the  vicinity  of  the  transitions  seems 
reasonable  based  on  a  study  of  the  surge  profiles  shown  in  Figures  20, 
21,  and  22.  In  addition,  prudent  design  would  suggest  some  minimum 
bottom  clearance  under  the  tow.  If  this  amount  is  set  at  0.5  ft,  which 
is  the  bottom  clearance  for  the  present  canal,  then  the  last  column  in 
Table  11  can  be  used  as  a  guide  in  alternate  canal  size  selection.  The 
table  shows  that  a  depth  of  15  ft  would  be  unnecessary  and  uneconomic. 
The  following  canal  sizes  give  approximately  equal  clearances  under  the 
tow. 


Approximate  Canal  Area 

Depth,  ft 

Width,  ft 

Increase,  % 

12 

275 

95* 

13 

200 

62* 

ll* 

Base 

39 

*  Interpolated  from  Table  5. 

62.  If  only  dredging  costs  are  considered,  the  most  economic 
solution  may  be  to  deepen  the  present  canal  to  ll*  ft.  However,  super¬ 
squat  effects  would  be  increased  from  the  present  canal  conditions  and 
the  increase  in  positive  surge  would  also  increase  the  tow  resistance 
in  the  transition.  The  increased  surge  conditions  may  cause  ar.  adverse 
reaction  from  navigation  interests  and  is  thus  not  recommended.  The  13- 
ft-deep  by  200-ft-wide  canal  size  would  not  increase  surge  magnitudes 
above  existing  conditions.  If  costs  can  be  justified  by  navigation 
traffic,  the  full  300-ft-wide  by  12-ft-aeep  canal  would  significantly 
improve  navigation  conditions.  The  depth  decrease  from  13  to  12  ft 
could  be  achieved  by  altering  maintenance  dredging  practice  if  the 
13-  x  200- ft  section  were  adopted  as  an  interim  improvement.  The 
recommended  canal  size  increases  would  have  the  following  impact  on  the 
main  parameters  involved  in  the  choice  of  canal  size. 
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Canal  Parameter 

Existing 

Unimproved 

Canal 

Base  Depth 
by 

3ase  Width 

Possible 

Intermediate 

Canal 

Improvement 
13-ft  Depth 
by 

200- ft  Width 

Eventual 

Canal 

as  Designed 
12- ft  Depth 
by 

300-ft  Width 

Canal  area,  sq  ft 

2000 

320C 

1»200 

Canal  to  tow  area  ratio 

2.2 

3.6 

h.l 

Approx  excav  vol,  M  cu  yd 

0.00 

0.20 

0.35 

Max  surge  amplitude,  ft* 

1-5 

1.3 

1.1 

Maximum  tow  speed,  mph 

2.h 

*».5 

5.3 

Steady  squat  at  max  speed,  ft 

0.7 

1.0 

1.0 

Supersquat ,  ft 

2.1 

2:i 

0.0 

Net  design  clearance,  ft** 

Grounding 

0.3 

1.0 

*  Obtained  from  Table  6. 

**  Obtained  from  Table  11,  assuming  1.0- ft  surge  amplitude  at  the  up¬ 
stream  canal  transition,  the  coincidence  of  tow  supersquat,  and  a 
0-5-ft  minimum  bottom  clearance. 


6 3.  The  following  conclusions  and  recommendations  are  based  on 
the  results  of  the  numerical  model  computations  presented  in  this  report. 
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a.  An  analysis  of  the  surge  computations  shows  that  canal 
widening  provides  greater  surge  reduction  than  canal 
deepening  for  a  given  increase  in  cross-sectional  area. 

The  magnitudes  of  the  reductions,  however,  are  modest, 
being  less  than  30  pcrc — '.t  decrease  of  the  present  surge 
with  100  to  150  percent  canal  area  increase. 

b.  Canal  area  increases  do  decrease  surge  magnitudes,  but 
the  largest  area  increase  to  more  than  twice  the  present 
area  cannot  bring  the  surge  below  about  1  ft. 

£.  Steady  squat  analysis  indicates  that  a  12-ft  depth  is 

the  minimum  necessary  to  eliminate  grounding.  Canal  width 
increases  without  deepening  would  be  detrimental  and  cause 
a  greater  tendency  for  tow  grounding. 

d_.  The  squat  analysis  shows  more  increase  in  limiting  tow 
speud  from  canal  deepening  than  widening. 

£.  Steady  squat  values  are  judged  to  be  acceptable  (less  than 
1  ft)  for  a  12-ft  canal  depth. 

f.  If  feasible,  speed  control  is  one  method  to  reduce  squat. 
However,  the  penalty  in  time  of  transit  is  not  attractive 
and  the  speed  limit  probably  could  not  be  enforced. 

£.  Increasing  the  canal  depth  without  width  increase  would 
increase  the  amount  of  supersquat  at  the  transitions. 
However,  the  net  underkeel  clearance  increase  is  greater 
than  the  supersquat  increase. 
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h.  A  200-ft-wide  by  12-ft-deep  channel  would  provide  about 
equal  relative  supersquat  to  the  present  canal  but  could 
still  result  in  grounding  under  adverse  surge  and  pool 
conditions.  A  12-  by  300-ft  canal  reduced  the  magnitude 
of  the  present  supersquat  to  steady  ahead  squat  values 
and  would  eliminate  the  possibility  of  grounding. 

i_.  Optimum  canal  to  tow  clearance  can  be  achieved  with  either 
of  the  following  minimum  canal  configurations: 


Dept1.,  ft  Width,  ft 


J  V 

i  l.  - 

12 

275 

4  4 

f 

1 

13 

200 

» 

.*> 

7 

J 

li< 

Base 

j_.  A  200-ft-wide  cana.'  with  a  minimum  13-ft 
511.0  to  reduce  grounding  problems  to  an 

depth  below  el 
acceptable  level 

•  i 
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would  provide  adequate  safety  with  minimum  dredging  cost. 

k.  It  is  recommended  that  canal  size  he  based  not  only  on  ex¬ 
cavation  cost  but  also  that  consideration  be  given  by  the 
Tulsa  District  to  maintenance  dredging  cost,  mobilization 
cost  for  future  required  widening,  and  benefits  from  greater 
tow  speed. 

1_.  If  navigation  traffic  can  justify  the  additional  cost, 
a  12-ft-deep  by  300-ft-wide  canal  should  be  dredged. 

This  would  significantly  improve  limiting  tow  speeds, 
probably  reduce  transit  times  through  the  canal,  and 
further  improve  navigation  conditions. 

m.  It  is  important  that  the  present  canal  not  be  widened 
without  deepening  as  this  would  aggravate  the  grounding 
problem. 
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Table  1 

Main  Cross-Sectional  Elements 
(at  Pool  El  311.0) 


Cross 

Section* 

Conveyance 
Top  Width 
ft 

Left  Bank 
to  0 

Nav  Chan 
ft 

Max 

Depth 

ft 

Storage 

Top 

Width 

ft 

Min 

Clearance  Under 
8-1/2-ft  Tow** 
ft 

1 

569 

178 

3*4.6 

1*4.5 

2 

522 

177 

28.9 

16. 5 

3 

511 

17*4 

2*4.2 

1*4.5 

b 

390 

121 

12.0 

1.8 

5 

370 

118 

12.3 

0.5 

6 

236 

122 

12.8 

2.0 

i 

232 

120 

12.0 

0.7 

8 

236 

120 

11.9 

0.6 

9 

380 

125 

11.3 

0.8 

10 

325+ 

55+ 

2*4.1 

12.5 

10-1 

2*43 

55 

2*4.1 

10-2 

177 

55 

2*4.1 

10-3 

110 

55 

25.0 

10— *4 

110 

55 

25.0 

RB-1 

UlU 

— 

27.8 

RB-2 

352 

— 

30.7 

LB-1 

3.8 

587 

LB-2 

3. *4 

11*4 

LB-3 

*4.  *4 

186 

LB- *4 

6.3 

*422 

Beat  ramp 

9.8 

355 

*  See  Figure  U  for  location  of  cross  sections. 

**  Tow  static  draft  of  8-1/2  ft;  105-ft  beam  tow  on  £  navigation  channel 
t  Section  to  face  left  lock  wall.  Total  section  width  approximately 
*468  ft. 


Table  2 

Summary  Table  of  Schematization 


Total 


Branch 

No. 

Length 

ft 

Grid  (Ax) 
ft 

No.  of 
Grids 

Name 

1 

2,985 

186.5625 

16 

River 

2 

1,257 

157.125 

8 

Spillway 

3 

l,Ul6 

177.00 

8 

Canal  entrance 

1» 

900 

75.00 

12 

1:2  transition 

5 

600 

60.00 

10 

Tow 

6 

2,97fc 

185.875 

16 

Canal 

7 

1,500 

150.00 

10 

1:10  transition 

8 

1,080 

180.00 

6 

Wide  canal 

9 

521 

130.25 

Lock  intake 

TOTAL 

13,233 

90 

_ Motes _ 

Upstream  boundary 
condition 

Upstream  boundary 
condition 

Tee  .junction 


Downstream  boundary 
condition 


Table  3 

Location  of  Model  Elements 

Station 

Cross 

ft 

Section 

Gage  Storage  Junction  Transition 

Notes 

River /Canal/Lock 


2227  +  90 

End  curve  downstream  bridge 

2229  +  10 

1 

1 

Upstream  model  limit 

2238  +  55 

7 

Downstream  old  Hwy  69  Bridge 

2239  +  35 

2 

22U1)  +  6!« 

1 

Beg  200'  It  bank 

22L6  +  61) 

1 

End  storage 

221)8  +  6!) 

3 

2258  +  80 

2 

Beg  1)50*  It  bank 

2258  +  95 

2 

?1) 1  depth 

Spillway  branch  No.  2 

to  12’ 

Rt  bank 

2261  +  83 

Beg  canal  -  300'  bottom  width 

2263  +  30 

2 

End  storage 

2265  +  20 

1) 

2269  +  1)1 

5 

300’  bottom  width  canal 

2270  +  11= 

Start  13°  curve  rt 

221)6  +  08 

221)8  +  98 

5 

221)9  +  08 

1) 

1:2 

Beg  300'  to  150'  bottom  width 

2252  +  08 

1:2 

End  300'  transition 

2252  +  3*» 

6 

2258  +  08 

5 

Canal/Tow 

Stern  tow  position 

2261)  +  08 

6 

Tow/Canal 

Bow  600 '  tow 

2268  +  83 

7 

2269  +  88 

End  13°  curve 

2273  +  33 

1) 

150'  bottom  width  canal 

2292  +  90 

8 

2293  +  82 

7 

1:10 

Beg  150'  to  300'  bottom  width 

2298  +  10 

3 

In  transition 

2308  +  82 

8 

1:10 

End  1500'  transition 

2309  +  05 

9 

2319  +  62 

9 

Canal/Lock 

End  300 '  bottom  width  canal 

Intake 

Beg  521 '  It  and  rt  eddy  zones 

2320  +  82 

Beg  It  guide  wall 

2320  +  92 

10-1 

>*,5 

Lt  and  rt  eddy  zones 

2372  +  00 

10 

2322  +  23 

10-2 

2322  +  71 

Beg  rt  guide  wall 

(Continued) 

Table  3  (Concluded) 


Table  5 

Summary  of  Pertinent  Canal  Cross-Sectional  Data 


No, 

.  8  for 

Pool  el  511.3) 

Nominal 

Dimensions* 

Area 

A 

sq  ft 

Top 

VgD 

ft/sec 

Relative 

Depth** 

ft 

Bottom 

Width 

ft 

Width 

W 

ft 

D  = 
A/W 
ft 

Area 

Increaset 

% 

1  “(an/a1) 

% 

Base 

Base 

1984.2 

238.4 

8. 32 

16.37 

0.00 

0.00 

Base 

200 

2412.7 

273.4 

8.82 

16.86 

17.76 

7.98 

Base 

250 

2880.7 

323.4 

8.91 

16.94 

45.18 

15.59 

Base 

300 

3347.7 

373.4 

8.97 

16.99 

68.72 

21.57 

12 

Base 

2342.0 

235-1 

9.96 

17-91 

18.03 

3.73 

12 

200 

2948.5 

280.1 

10.53 

18. 4l 

48.60 

13.01 

12 

250 

3563.5 

330.1 

10.80 

18.64 

79.59 

20.37 

12 

300 

4178.5 

380.1 

10.99 

18.81 

no.  59 

26.13 

14 

Base 

2760.2 

238.4 

11.58 

19.31 

39.11 

7.92 

Ik 

200 

3476.2 

289.4 

12.01 

19.67 

75.19 

17-19 

l4 

250 

4191.2 

339.4 

12.35 

19.94 

111.23 

24.06 

15 

Base 

2999-4 

242.4 

12.37 

19-96 

51.16 

10.19 

15 

200 

3764.4 

292.4 

12.87 

.  >.36 

89.72 

19-04 

15 

250 

4529.4 

342.4 

13.23 

20.64 

120.27 

25.69 

15 

300 

5294.4 

392.4 

13.49 

20.84 

166.83 

30.92 

*  Base  (existing)  depth  is  9  ft  and  base  width  is  150  ft. 

**  Nominal  depth  at  normal  upper  pool  el  511.0. 
t  Computed  for  each  canal  configuration  compared  with  existing  canal 
area  of  1984.2  sq  ft. 


Table  6 

Summary  of  Surge  Calculations  at  Various  Canal  Conditions 


4g 


Canal 

Test 

Cond 

Min  £1 
§  6  Min 
ft  ms  3 

Surge 

Amp 

ft 

Surge  Amp 
Reduction 
ft 

Rel  Amp 
Reduction 
% 

Min  El  § 
Sta  77/78 
ft  msl 

Surge 

Amp 

ft’ 

Surge  Amp 
Reduction 
ft 

Rel  Amp 
Reduction 
% 

;  1 

Base  x 

Base 

509-96 

1.3b 

0.00 

0.00 

509.85 

l.b5 

0.00 

0.00 

1 

Base  * 

200 

510.06 

1.2b 

0.10 

7.b  6 

509-98 

1.32 

0.13 

8.97 

3 

Base  x 

250 

510.17 

1.13 

0.21 

15.67 

510.09 

1.21 

0.2b 

3  5.55 

Base  x 

300 

510.25 

1.05 

0.29 

21.6b 

510.19 

1.11 

0.3b 

23. b5 

^4 

12  x 

Base 

509-93 

1-37 

-0.03 

-2.2b 

509-86 

l.bb 

0.01 

0.7 

12  v 

200 

510.06 

1.2b 

0.10 

7. bo 

510.00 

1.30 

0.,S 

10.3b 

1  3 

i  4g 

12  » 

25c 

510.18 

1.12 

0.22 

I6.b2 

510.12 

1.18 

0.27 

18.62 

{  2 

12  x 

300 

510.25 

1.05 

0.29 

21.6b 

510.21 

1.09 

0.36 

2b.  83 

M 

3 

lb  x 

Base 

510.0b 

1.26 

0.08 

5.97 

509-97 

1.33 

0.12 

8.28 

1 

lb  x 

200 

510.1b 

1.16 

0.18 

13.b3 

510.10 

1.20 

0.25 

17.2b 

-z* 

% 

lb  x 

250 

510.23 

1.07 

0.27 

20.15 

510.20 

1.10 

0.35 

21.1b 

1 

15  x 

Base 

510.07 

1.23 

0.11 

8.21 

,10.01 

1.29 

0.16 

11.03 

f  M 

15  x 

200 

510.19 

2 . 11 

0.2.J 

17.16 

510. lb 

3.16 

0.29 

20.00 

15  x 

250 

510.28 

1.02 

0.32 

23.88 

510.29 

1.01 

O.bb 

30.3b 

'  ^ 

■ 

15  * 

300 

510.35 

0-95 

0.39 

29-10 

510.36 

0.9b 

0.51 

35.17 

* 

i  m 
!  M 

f 


?  ■  ' 

|  U4 

•4 

■  % 
«  1 4 

t  »■} 


Table  7 

Summary  of  Souat  Calculations  at  Various  Canal  Conditions 
(Geometry  at  Section  8  -  See  Table  5) 


Canal  Test 
Conditions 

Limiting 

Tow 

Speed 

Vt 

ft/s ec 

Relative 

L ini ting 
Tow  Speed 
Increase 

* 

Return 

Velocity 

ui 

ft/sec 

Squat 

at 

vt 

ft 

90S  X 
Vt 

ft/sec 

Squat 
at  0.9V t 

ft 

Relative 
Increase 
in  Squat 
at  0.9Vt 

S 

Constant 

Speed* 

Squat 

ft 

Relative 
Constant 
Speed  Squat 
Decrease 

S 

Base  x  Base 

3-925 

0.00 

6.2L5 

1.367 

3-533 

0.697 

00.00 

0.697 

00.00 

Base  x  200 

5-098 

29. 87 

6.217 

1.588 

8.588 

0.793 

13.82 

0.366 

87.55 

Base  x  250 

6.031 

53.63 

5.973 

1.673 

5-827 

0.822 

17.99 

0.288 

68.81 

Base  x  300 

6.780 

72.72 

5.729 

1.716 

6.102 

0.829 

19-03 

0.188 

73.01 

12  *  Base 

5-289 

33.73 

6.61*7 

1.770 

8-725 

0.889 

97  r,o 

0.381 

85.33 

12  x  200 

6.681 

70.20 

6.L51 

1.985 

6.013 

0-973 

39-63 

0.231 

66.83 

12  x  250 

7  759 

99.67 

6.160 

2.078 

6.983 

0.995 

82.81 

0.166 

76.22 

12  x  30c 

8.618 

119-51* 

5.885 

2.113 

7-756 

O.998 

82.70 

c.129 

81.87 

lU  x  Base 

6.311 

68-93 

6.8SO 

2.156 

5-968 

1.065 

52.78 

0.259 

62.85 

lh  X  200 

8.052 

105.13 

6.551 

2.305 

7.287 

1.109 

59-18 

0.171 

75.88 

1*.  x  250 

9.189 

133.07 

6.231 

2-373 

8.238 

1.116 

60.22 

0.127 

81.71 

15  *  Base 

7.31*3 

87.06 

6.959 

2.339 

6.609 

1.185 

68.  ?6 

0.220 

68.89 

15  *  200 

6.773 

123-50 

6.605 

2.877 

7.896 

l.!-1! 

60.87 

0.150 

78.58 

a  5  *  250 

9-870 

151-1*1* 

6.270 

2-532 

5.883 

1.179 

69-27 

0.113 

83.78 

15  x  300 

10.71*3 

173.67 

5.969 

2.585 

9.668 

I.I61 

66.62 

0.C912 

86.92 

*  Constant  speed  used  is  equal  to  0.9vj  Cor  the  existing  (o'  *  150')  canal  which  is  2. hi  sph. 


'  VMM  I  *  4«,l  in  fHU  In . I  III  Mlttl  I.  tigs . .  MM  . . .  ,m  ny,.,,,  ,|,  . . . . 


Tabic  8 

Sussarg  of  Supersquat  Calculations  at  1:2  Transition 
for  12- ft  Depth  and  Various  Widths 


I  !? 


Table  10 

of  Supersquat  Calculations  at  1:2  Transition 
For  lk-ft  Depths  and  Various  Widths 


lfr*  x  Base 


Ik'  x  200* 
Section  Numbers 


Ik*  x  250* 


Variable 


U955.0  2772.1  1*955.0  3k8Y-l  **955-0  **202.1 


w 

396.5  21*2.1 

396.5 

292.1 

396.5  3k2.1 

n  =  A/a* 

5.552  3.106 

5-552 

3.907 

5.552 

k.708 

D  =  A/W 

12.1*97  11.1*50 

12.k97 

11.938 

12.k97 

12.283 

Vg5 

20.0o0  19.201 

20.060 

19.606 

20.060 

19.887 

UJt 

5.892  6.81*1* 

5-892 

6.525 

5-892 

6.208 

Z£ 

2.37**  2.13k 

2.37k 

2.291 

2.37k 

2.360 

Vo 

10.028  6.619 

10. 028 

8.0Uk 

10.028 

9-138 

0-9V, 

9-025  5-957 

9-025 

7.2k0 

9.025 

8.22k 

Z0-9V£ 

1.093  1-053 

1.093 

1.102 

1.093 

1.110 

uo.9V 

£ 

3.297  k.208 

3-297 

3.869 

3-297 

3-571 

0.9V .  ^ 
*© 

© 

> 

1 

2.ko6 

0.981 

-0.113 

n»* 

l.k35 

0.597 

0.00 

Ah  =  n  + 

l£©f 

3.569 

2.888 

1.861 

ah 

»: 

"  "  Ahl 

-33.36 

-7-93 

30.  k7 

*  a  =  design  tow  cross-sectional  area  =8.5  ft  x  105  ft  =  892.5  sq  ft. 
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If  n  <  0  ,  Ah  =  Vt  =  V9Vj 
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Table  11 

Summary  of  Tow  Clearance  Calculations 
at  Various  Canal  Conditions 
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Canal  1 

Size 

Design  Clear¬ 
ance  Under 
8-1/2-ft  Tow 
ft 

Supers quat 
Ah 
ft 

Supersquat 

Clearance 

ft 

Net  Design 
Clearance* 
ft 

Base  x 

Base 

0.5 

2.7 

Grounding 

Grounding 

Note: 

Increasing  canal 

width  without  depth  increase  does  not 

prevent  grounding  at  steady  or 

supersquat 

tow  speed.’. . 

12’  x 

Base 

3.5 

3.2 

0.3 

Possible  grounding 

12’  x 

200' 

3.5 

2.6 

0.9 

Possible  grounding 

12*  x 

250' 

3.5 

2.1 

1.1* 

Marginal 

12'  x 

300' 

3.5 

1.0 

2.5 

1.0 

13'  x 

Base 

1*.5 

3.1*** 

1.1 

Marginal 

13'  x 

200' 

U.5 

2.7** 

1.8 

0.3 

1*+'  X 

Base 

5.5 

3.6 

1.9 

0.1* 

lV  X 

200' 

5-5 

2.9 

2.6 

1.1 

lh<  x 

250' 

5-5 

1.9 

3. 6 

2.1 

15'  x 

Base 

6.5 

3.8 

2.7 

1.2 

15'  x 

200' 

6.5 

3.0 

3.5 

2.0 

15'  x 

250' 

6.5 

1.8 

1*.  7 

3-2 

15:  x 

300' 

6.5 

1.2 

5-3 

3.8 

M 


o 

t 


J 


I 

ft '  J 
'  <? 


1 !  I 
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*  Supersquat  clearance  minus  lock  surge  allowance  (1.0  ft)  and 
minimum  bottom  clearance  (0.5  ft). 

**  Linearly  interpolated  supersquat. 
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APPENDIX  A:  NOTATION 


Tow  cross-sectional  or  frontal  area,  sq  ft 

Surge  amplitude  for  the  present  unmodified  canal,  ft 

Surge  amplitude  for  any  given  variable  canal  geometry,  ft 

Cross-sectional  flow  area  of  channel,  sq  ft 

Celerity  of  surge  in  a  channel  in  ft/sec;  for  small  amplitude 
surges ,  C  =  VgD 

Mean  depth  of  channel,  ft,  5  A/W 

Af/Wi  present  canal  mean  depth,  ft 

Variable  geometry  canal  mean  depth,  ft 

2 

Acceleration  due  to  gravity,  ft/sec 

Water-su. face  elevation  above  mean  sea  level,  ft 

Manning's  resistance  coefficient  for  surge  computations; 
ratio  of  channel  area  to  tow  area  for  tow  squat  computations 

Lateral  inflow  per  unit  distance  along  the  channel  per  unit 
time,  sq  ft/sec 

Total  channel  discharge,  cu  ft/sec,  normally  Q  =  AV 
Hydraulic  radius  of  channel,  ft 
Return  velocity,  ft/sec 

Initial  estimate  of  return  velocity  for  steady  squat  compu¬ 
tations,  ft/sec 

Return  velocity  at  limiting  tow  speed,  ft/sec 
k+l  estimate  by  iteration 
Mean  flow  velocity  of  channel,  ft/sec 
Schi.jf  limiting  tow  speed,  ft /sec 

Schijf  limiting  velocity  at  given  channel  cross  section  (?)  , 
ft/sec 


A1 


Tow  speed,  ft/sec 


Tow  speed  at  given  cross  section  (!)  for  calculation  tow 
supersquat,  ft/sec 

Top  width  of  water  surface  of  channel,  ft 

Present  canal  top  width,  ft 

Variable  geometry  canal  top  width,  ft 


I  I 'I 


Vertical  tow  squat,  ft 

Initial  estimate  of  tow  squat  for  steady  squat  computations, 
ft 

k+1  estimate  by  iteration 

Tow  squat  ab  limiting  tow  speed,  ft 

Magnitude  of  tow  supersquat  at  an  abrupt  channel  transition, 
ft;  =  z£  +  n 

Time-step  for  lock  surge  model  computations,  sec 

Space-step  for  lock  surge  model  computations,  ft 

Positive  surge  magnitude  induced  by  a  tow  driving  through 
an  abrupt  channel  transition,  ft 

Rate  of  change  with  respect  to  time,  per  time 


Rate  of  change  with  respect  to  distance,  per  ft 

90  percent  of  limiting  tow  speed  at  cross  section  (l)  , 
which  is  the  normal  upper  possible  speed  of  self-propelled 
tows,  ft /sec 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Huval,  Carl  J 

Lock  approach  canal  surge  and  tow  squat  at  Lock  and 
Dam  17,  Arkansas  River  project;  mathematical  model 
investigation  /  by  Carl  J.  Huval.  Vicksburg,  Miss.  : 

U.  S.  Waterways  Experiment  Station  ;  Springfield,  Va.  : 
available  from  National  Technical  Information  Service, 

198° . 

53,  [12],  2  p.  :  ill.  ;  27  cm.  (Technical  report  -  U.  S. 
Army  Engineer  Waterways  Experiment  Station  ;  HL-80-17) 
Prepared  for  U.  S.  Army  Engineer  District,  Tulsa, 

Tulsa,  Oklahoma. 

References:  p.  53. 

1.  Arkansas  River.  2.  Canalization.  3.  Lock  filling  and 
emptying  systems.  L.  Locks  (Waterways).  5.  Mathematical 
models.  6.  Navigation  channels.  7.  Squat  of  vessels. 

8.  Surges.  I.  United  States.  Army.  Corps  of  Engineers. 
Tulsa  District.  II.  Series:  United  States.  Waterways 
Experiment  Station,  Vicksburg,  Miss.  Technical  report  ; 
HL-80-17. 
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